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Abstract 

A Gaussian multiple-input single-output (MISO) wiretap channel model is considered, where there 
exists a transmitter equipped with multiple antennas, a legitimate receiver and an eavesdropper each 
equipped with a single antenna. We study the problem of finding the optimal input covariance that 
achieves ergodic secrecy capacity subject to a power constraint where only statistical information about 
the eavesdropper channel is available at the transmitter. This is a non-convex optimization problem that is 
in general difficult to solve. Existing results address the case in which the eavesdropper or/and legitimate 
channels have independent and identically distributed Gaussian entries with zero-mean and unit-variance, 

f*"*- , i.e., the channels have trivial covariances. This paper addresses the general case where eavesdropper and 

CO ■ 

■ legitimate channels have nontrivial covariances. A set of equations describing the optimal input covariance 

matrix are proposed along with an algorithm to obtain the solution. Based on this framework, we show 

o 

that when full information on the legitimate channel is available to the transmitter, the optimal input 
covariance has always rank one. We also show that when only statistical information on the legitimate 
channel is available to the transmitter, the legitimate channel has some general non-trivial covariance, and 

X! ■ 

the eavesdropper channel has trivial covariance, the optimal input covariance has the same eigenvectors 



C3 



as the legitimate channel covariance. Numerical results are presented to illustrate the algorithm. 

Index Terms 

Ergodic secrecy capacity, MISO wiretap channel, beamforming. 
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I. Introduction 

Wireless physical (PHY) layer based security from a information-theoretic point of view has received 
considerable attention recently |1]. Such approaches exploit the physical characteristics of the wireless 
channel to enhance the security of communication systems. The wiretap channel, first introduced and 
studied by Wyner |2], is the most basic physical layer model that captures the problem of communication 
security. Wyner showed that when an eavesdropper's channel is a degraded version of the legitimate 
channel, the source and destination can achieve a positive information rate (secrecy rate). The maximal 
secrecy rate from the source to the destination is defined as the secrecy capacity; for the degraded wiretap 
channel the secrecy capacity is given as the largest between zero and the difference between the capacity 
at the legitimate receiver and the capacity at the eavesdropper. The Gaussian wiretap channel, in which the 
outputs at the legitimate receiver and at the eavesdropper are corrupted by additive white Gaussian noise 
(AWGN), was studied in 0. Along the same lines, the secrecy capacity of a deterministic Gaussian MIMO 
wiretap channel has been studied recently in E1|-|S). In (9J, the achievable rate in Gaussian MISO channels 
was studied. In that context, the channel state information (CSI) of the legitimate channel was assumed to 
be available, but only statistical information about the eavesdropper channel was assumed to be available 
at the transmitter. In (9] it was shown that when the eavesdropper channel is a vector of independent and 
identically distributed (i.i.d.) zero-mean complex circularly symmetric Gaussian random variables, i.e., 
the channel has a trivial covariance matrix, the optimal communication strategy is beamforming, and that 
the beamforming direction depends on the CSI of the legitimate channel. In iflOl , the authors derived the 
ergodic secrecy capacity of a Gaussian MIMO wiretap channel where only statistical information about 
the legitimate and eavesdropper channels are available at the transmitter. It was shown that a circularly 
symmetric Gaussian input is optimal. It was also shown in the same paper that when the eavesdropper 
and legitimate channels have i.i.d. Gaussian entries with zero-mean and unit-variance (trivial covariance), 
a circularly symmetric Gaussian input with diagonal covariance is optimal. 

In this paper, we consider a Gaussian multiple-input single-output (MISO) wiretap channel and assume 
that only statistical information about the the eavesdropper channel is available at the transmitter. Re- 
garding the legitimate channel, we consider two scenarios: a) only statistical information of the legitimate 
channel is available at the transmitter; b) full CSI on the legitimate channel is available at the transmitter. 
We extend the result of [9] and [ 10] proposed for the case of multiple-input single-output (MISO) wiretap 
channel with trivial channel covariances to the case of nontrivial covariances. The non-trivial channel 
covariance matrix corresponds to the case where there exists statistical correlation between the channel 
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coefficients of different transmit-receive antenna pairs. Such cases arise when the transmit and receive 
antennas are closely spaced relative to the signal wavelength. We address the problem of finding the 
optimal input covariance that achieves ergodic secrecy capacity subject to a power constraint. This leads 
to a non-convex optimization problem. The contributions of this paper are the following: 

• We derive a set of equations for the optimal input covariance matrix, and propose an algorithm to 
obtain the solution (please refer to Theorem 1 of Section IV). 

• We show that when the legitimate channel is completely known at the transmitter, in addition to the 
conditions of Theorem [T] the following hold: 1) the optimal input covariance matrix has rank one; 
2) the ergodic secrecy rate is increasing with the signal-to-noise ratio (SNR). 

• We show that when only statistical information on the legitimate channel is available to the transmit- 
ter, the legitimate channel has some general non-trivial covariance, and the eavesdropper channel has 
trivial covariance, the optimal input covariance has the same eigenvectors as the legitimate channel 
covariance. 

• We show that under high SNR, the optimal input covariance has rank one. 

The remainder of this paper is organized as follows. The mathematical model is introduced in Spl] 
In fllTTl we give the explicit expression of ergodic secrecy rate, and in SJTV] we derive the condition for 
optimal input covariance. In fJVJ we analyze the dependence of ergodic secrecy rate on the SNR, and in 
WH we study the ergodic secrecy rate under high SNR. In Will an algorithm is proposed to search for 
the solution. Numerical results are presented in Willi to illustrate the proposed algorithm. Finally, E jlXl 
gives a brief conclusion. Several proofs appear in an Appendix. 

A. Notation 

Upper case and lower case bold symbols denote matrices and vectors, respectively. Superscripts *, 
T and f denote respectively conjugate, transposition and conjugate transposition. det(A) and Tr(A) 
denote the determinant and trace of matrix A, respectively. A max (A) and A m i n (A) denote the largest and 
smallest eigenvalues of A, respectively. A y means that A is Hermitian positive semi-definite, and 
A y means that A is Hermitian positive definite, diag(a) denotes a diagonal matrix with the elements 
of the vector a along its diagonal. ||a|| denotes Euclidean norm of vector a. I n denotes the identity matrix 
of order n (the subscript is dropped when the dimension is obvious). E{-} denotes expectation operator. 
In this paper, log(-) denotes base-e logarithm where e = 2.71828 
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Fig. 1. System model. 

II. System Model and Problem Statement 

Consider a Gaussian MISO wiretap channel shown in Fig. [T] where the transmitter is equipped with 
tit antennas, while the legitimate receiver and an eavesdropper each have a single antenna. The received 
signals at the legitimate receiver and the eavesdropper are respectively given by 



where x is the rifX 1 transmitted signal vector with zero mean and ny x nj- covariance matrix R x . y 0, 
i.e., x ~ CAf(0, R x ); h.R, Iie are respectively channel vectors between the transmitter and legitimate 
receiver, and between the transmitter and eavesdropper; v n ~ CAf(0, <j%), ve ~ CM(0, a^) are the noises 
at the legitimate receiver and the eavesdropper, respectively. We can represent H x in terms of the average 
signal energy E s and normalized signal covariance matrix Q, so that H x = E S Q and Tr(Q) = 1. The 
signal-to-noise ratio (SNR) is defined as p = E s / 'a^. 

We assume that full CSI is available at both the legitimate receiver and the eavesdropper, and only 
the statistical information on the eavesdropper channel is available at the transmitter. We consider two 
cases, depending on the type of information available at the transmitter on the legitimate channel: 
a) Only statistical information on the legitimate channel is available at the transmitter, i.e., the trans- 
mitter knows the distributions of h# and given by ~ CAf(0, he ~ CJ\f(0, T,e) with 
covariances y 0, and "Ee y 0, respectively. The ergodic secrecy capacity of the Gaussian MISO 
wiretap system Q equals [10] 



VR = h i? x + V R 
y E = h^x + v E 



(1) 



(2) 



a 



A 




(3) 



s 
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where C S (Q) is the ergodic secrecy rate given by 



C s (Q)=E hfl {log(l + phkQh fl )} 



E hE {log(l + phT,Qh E )}. 



(4) 



b) Full CSI on the legitimate channel is available at the transmitter. The ergodic secrecy rate is given 



The transmitter optimization problem is to find the optimal input covariance matrix Q to maximize 
C S (Q) for cases a) and b). We denote the feasible set as SI = {Q[Q h 0, Tr(Q) = 1} which is a convex 
set. 

The problem is of interest when a positive secrecy rate can be achieved, i.e., C S (Q) > for some Q. 
The conditions to ensure a positive ergodic capacity are provided in the following lemmas. 

Lemma 1: For hf( ~ CAf(0, S_r), the sufficient and necessary condition under which C S (Q) > for 
some Q is that 5]^ — He is non negative semi-definite. 
The proof is given in Appendix |A] 

Lemma 2: When h/j is completely known at the transmitter, a sufficient condition under which 
C S (Q) > for some Q is that hRh^ — He is non negative semi-definite. 
The proof is given in Appendix iBl 



The calculation of the ergodic secrecy rate involves calculation of terms like E z {log(l + /?z^Qz)} with 
z ~ CAA(0,R). To this end, following the analysis of |[TT1 Eq. (64)], we give the following lemma. The 
proof is given in Appendix ICl 

Lemma 3: Let R 1 / 2 QR 1 / 2 have a eigen-decomposition UiDiUj where Di = diag((fi, • • • , dm , 0, • • • , 0), 
and d\ > ■ ■ ■ > dM > are the M non-zero eigenvalues. For z ~ CAf(0, R), it holds: 



by El 



C a (Q) = log(l + ph^Qhfl) - E hE {log(l + P h E Qh E )}. 



(5) 



III. Calculation of Ergodic Secrecy Rate 




M 



(6) 



where Fi(x) 



e l l x ~Eii(l/x) with Ei(x) = ^j-dt being the exponential integral. 



Based on (J6]), we can calculate C S (Q) by simply letting R = or R = S^. 
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IV. Conditions for Optimal Input Covariance 

Next we obtain the necessary conditions for the optimal Q by using the Karush-Kuhn-Tucker (KKT) 
conditions. Let us construct the cost function 

L(Q, 6, *) = C S (Q) - 0(Tr(Q) - 1) + Tr(*Q) (7) 

where 6 is the Lagrange multiplier associated with the constraint Tr(Q) = 1, and ^Sf is the Lagrange 
multiplier associated with the constraint Q >z 0. The KKT conditions enable us to write EDI 



© - ei nT + * = o, (8) 

* h 0, Tr(*Q) = 0, Q h 0, Tr(Q) = 1, (9) 
where = (^p) T . By using the fact = (hht) T , we have: for case a) 

0=E h / Ph f« \-eJ Ph t* ] (10) 



lR 1 1 + ph^Qh R J E 1 1 + ph E Qh E 
and for case b) 

0=^# — eJ ph f° y (id 

From the KKT conditions ([8]) and (O, we obtain the equivalent (but without containing the Lagrange 
multipliers) conditions for optimal Q consisting of a set of equations given in the following theorem. 
Please see Appendix O for details. 

Theorem 1: The optimal Q >z satisfies 

Q0 = 0Q = Tr(Q0)Q (12) 

A max (0) = Tr(Q0). (13) 

The above conditions imply that for the optimal Q, Q0 is a scaled version of Q. Any Q satisfying the 
conditions of Theorem Q] is called as KKT solution. In Will we propose an algorithm to search for the 
KKT solution. is an important variable for the transmitter optimization problem. For the calculation 
of 0, we give the following lemma. The proof is given in Appendix El 

Lemma 4: For z ~ CM(0, R), it holds 

E 4t^TFt) = pR^UiYUtR 1 ^ (14) 
[ 1 + pzTQz J 
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where Y is a diagonal matrix with (k, k)th entry given by 



M 



Y kk = Yl 



n5»(i-*/*) 

M 



+ _ M ; K K : — ,k<M (i5) 



Y kk =Y j ,k>M. (16) 

i= i pdjUi^ji 1 -di/dj) 

with F2(x) = - — ^2-e 1//a: Ei(l/a;), and F\(x), E\(x), Ui, di, M are defined in Lemma\3\ 
Based on (fl4l) . we can calculate by simply letting R = or R = Sg. 

In the following, we show that for some special cases, more information about Q than that of Theorem 
1 can be obtained. 



A. h/j is completely known at the transmitter 

We will prove that if C S (Q) > for some Q, then the optimal Q always has rank one, i.e., beamforming 
is optimal. We put the proof in the second part of the subsection. In the first part of the subsection, we 
analyze how this result reduces our problem to a problem of one variable. 

Based on this result, we let Q = uu^ with ||u|| 2 = 1 and the problem is reduced to 

C S (Q) = logfl + phJjimthjO -E hE {log(l + p4uuth s )} (17) 

which, by using (l66l ) and (I68T ). can be rewritten as 

C S (Q) = log(l + puth^u) - Fi(putS^u). (18) 

Let u^h^h^u = z||h/j|| 2 . Then < z < 1. Note that F\(x) is an increasing function. Thus, for fixed z, 
u^eu should be minimized. Define 

4>(z) =min m^em (19) 

u 

s.t. u^h^hJjU = z||h#|| 2 , and ||u|| 2 = 1. 

Then, our problem is reduced to 

C,(z) = log(l + p||h R || 2 z) - F 1 (p<P(z)), 0<z<l. (20) 
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Since the problem of ( fl9l ) belongs to the class of quadratically constrained quadratic programming 
(QCQP) with two constraints, it can be exactly solved iTTTl . and is equivalent to its semidefinite program- 
ming (SDP) relaxation, i.e., 

<f>(z) =min TV(S £ X) (21) 

s.t. Tr(h R h^X) = ^ 1 1 1 1 2 , andTr(X) = 1, 

X >- 0. 



For any given z, the problem of (12Tb is an SDP and can be effectively solved via CVX software [25]. 

Lemma 5: The function <p{z) is a convex function. 
The proof is given in Appendix IB 

Since <fi{z) is a convex function, according to well-known properties of convex functions, we know 
that <ft(z) is continuous and Lipschitz continuous |[27l Corollary 2.3.1], and is differentiable at all but at 
most countably many points (left and right derivatives always exists) [27 , Theorem 2.3.4]. Further study 
on cf){z) and proposing more effective method for the optimization of C s (z) can be our future work. 

For the special case S# = al (hg has a trivial covariance), (fT8l) becomes 

C S (Q) = log(l + piA^u) - Ftipa). (22) 

Obviously, the optimal u = hfl/||hfl||, the optimal Q = h/jh^,/||h/j|| 2 and 

[maxC s (Q)] = log(l + p\\h R \\ 2 ) - F^pa). (23) 

This case was considered in |9] and the above result is consistent with that in |9ll . 
Remark: If the optimal u do not achieve C S (Q) > 0, then C S (Q) < for any Q. 

In the remainder of the subsection, we give the proof for that if C S (Q) > for some Q, then the 
optimal Q always has rank one. We first provide a lemma that will be helpful in the following. The proof 
is put in Appendix iGl 

Lemma 6: Let A be a positive definite matrix, a. be a vector. If aa^ — A has a positive eigenvalue, 
then it has all negative eigenvalues except for a positive eigenvalue. 

Via Lemma [6l we can show that, if C S (Q) > 0, then has all negative eigenvalues except for a 
positive eigenvalue. To see why this is the case, recall that 



@ = ^^ Ej Ph f* }. (24) 

1 + ph R Qh R 1 1 + ph E Qh E . 
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According to W- II (after Lemma 01, we know that, if C S (Q) > 0, then Tr(0Q) > 0. Thus, has at 
least a positive eigenvalue. According to (fT4b . we know that the second term in the right hand side of 
(1411) is positive definite. Note that the first term in the right hand side of (1411) has the form aa^. Thus, 
the desired result follows directly from Lemma [6] 

From Theorem [T] we know that the optimal Q and its associated are commutable. Thus, there exists 
a unitary matrix Uo that simultaneously diagonalizes Q and 0. Let Aq and A® be the corresponding 
diagonal matrices. From (TT2l) in Theorem [T] we know that 

A Q A@ = Tr(0Q)A Q (25) 

or equivalently, 

(A Q ) kk (A & ) kk = Tr(0Q)(A Q ) fefc , k = l,---,n T . (26) 

Since Tr(0Q) > 0, it follows from (|26]> that, if (A Q ) kk > 0, then (A & ) kk = Tr(0Q) > 0. However, 
A© has all negative diagonal entries except for a positive one. Thus, Aq has only one nonzero diagonal 
entry, i.e, the optimal Q has rank one. 

B. Only statistical information on \ir available at the transmitter and He = oil 

During this subsection, we assume that E# has simple spectrum (all eigenvalues are distinct), since mul- 
tiple eigenvalues are rare for generic Hermitian matrices ['19, §4]. Let have the eigen-decomposition 
£r = V^A/jV^ where A R = diagfai,-- - ,% T ), rji > 7] 2 > ■ ■ ■ > i] nT - We use Theorem Q] to 
show that the optimal Q has the same eigenvectors as Xr, i.e., vJjQVjj is diagonal, denoted by 
A = diag(£i, • • • , Cn T )- Using this result, our problem is reduced to 

C S (Q) = E h Jlog(l + ph^ A]( 2 A A^ 2 h w ) } 

-E hw {log(l + apht ) Ah 11) )} 

= E hu ,{log(l + p V" T rjid\h w .i\ 2 )} 

* J l=l 

. Xi\h w>i 2 )}. (27) 

The power constraint is Yl7=i & = ^ ^ or tne case nT = ^» < l27l ) becomes 

C.(Q) = E hm {log(l + pviCiKM 2 + PV2(1 - Ci)l^, 2 | 2 )} 

- E h „{log(l + apd\ Vil 2 + ap(l ~ Cl)l V2| 2 )}- (28) 

The constraint is < £i < 1. Similarly to S JTITI the expectations in (|27T ) and ([28]) can be expressed in 
explicit form. 
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In the remainder of the subsection, we give the proof. Let Q have the eigen-decomposition Q = VAV^ 
where A is diagonal. According to Appendix [0 we express as 

= Eh { p4\^4 \ _ Eh J ap ^ ) . (29) 



1 + phts^Q'S 1 ^ J " 1 1 + aphlQh w ) 

i i t 
Let SjjQSjj = A and let A have the eigen-decomposition A = V A A A V A where 

A A = diag(Ai Ji, A2J2, • • • , X k Jk), (30) 

Ai > A2 > • • • > Xk-i > Ak > are distinct eigenvalues, and J^'s are identity matrices. By using the 
fact that Uha, and h w have the identical distributions for any unitary matrix U, we express as 

= Y^YaY^Y^ 2 - VZ Vt (31) 



where 



■ ■ { Ph f W 1 (32) 
1 + ph^A^h 

Eu< aphwhl . \. (33) 



1 + aphJuAh. 

Similarly to Appendix El it can be shown that Yo and Zq are both diagonal. 

Observe that VZoVt and Q are commutable. With this, from Theorem [T] we know that and Q are 
commutable which enables us to get 

E^VaYqV^S^Q = QS^VaYqV^ 2 . (34) 

By inserting EJjQEJg = V A A A V A into d34l . we get 

S r VaYoAaV a = V A A A Y V A S R . (35) 

Since A A and Yo are both diagonal matrices, it holds that YoA A = A A Yo. With this, by inserting 
S fl = VrArV^ into OD, we get 

A R V A A Y V = V A A Y vlA R (36) 

where Vo = V^V A . From (l36l) and the assumption that all diagonal entries of Ar are distinct, we know 
that VoA^YoVj is a diagonal matrix |[22l Special matrices: diagonal]. On the other hand, it follows 
from £|.Q£! = V A A A v\ that Q = S^ 1/2 V a AaV^£^ 1/2 which, when combined with the fact that 
Sr = VijAflVL results in 



Q = V R A R 1/2 (V A A vl)A- 1/2 V j R . (37) 
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Next, we show that VoA^vJ is diagonal. Since VoA^YoVq is diagonal, there exists a Vi, which 
is the column-permuted version of Vo, such that ViA^YoVj = A^Yo, or equivalently, ViA^Yo = 
A A Y Vi. We aim to prove that ViA^ = A^Vi. According to ([30]>, (1321 ). ( |76l ), ( 1771 ), it is not difficult 
to show that A^Yo has the form of 

A A Y = diag(Ai Ji, \' 2 J 2 , ■■■ , \' k Jk) (38) 

where X[ > X' 2 > ■ ■ ■ > A' x _ x > 0, \' K > 0, and \' K ^ A^ for k = 1, • • ■ , K - 1. From @8]l and the fact 
that ViA^Yo = A^YoVi, we know that Vi has the form of Vi = diag(Ai, A2, ■ ■ ■ , Ax), where each 
Afe is the same size as the corresponding ll22l Special matrices: diagonal]. Thus, it is easy to verify 
that ViA^ = A4V1, or equivalently, ViA^V"} = A^. Therefore, since Vi is the column-permuted 
matrix of Vo, it follows that VoA^vJ is diagonal. With this, from (|37T ). we know that the optimal Q 
has the same eigenvectors as S/j. 

V. Dependence of C s (Q) on p 

In this section we investigate how the SNR, p, impacts the ergodic secrecy rate. 
1 ) Full CSI on h.R at the transmitter: We first provide a lemma that will be helpful in the following. 
The proof is given in Appendix [H] 

Lemma 7: For a positive constant x and a positive random variable Y, the following fact holds: 

logx > E(logF) => - <E ( — ) . (39) 



x \Y, 

Here, ==>■ means that the right side follows from the left side. 

By using Lemma |7J we can prove that, if C S (Q) > 0, then Tr(@Q) > 0. To see why this is the case, 
we let x = 1 + ph^QhR and Y = 1 + ph^Qhg which enables us to write 

C S (Q) = logx-E(logY) (40) 

Tr(0Q)=Efi^ -i. (41) 

The desired result follows from Lemma |7J 

Taking the derivative of C S (Q) with respect to p, we get: 

dC s (Q) _ h R Qh R f h E Qh E 

dp " 1 + ph R Qh R hE 1 1 + ph^Qh s 

= . (42) 

P 
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Based on the fact that, if C a (Q) > 0, then Tr(0Q) > 0, we get that if C S (Q) > 0, then 8C ^ Q) > 0. 
Thus, if C S (Q) > for some Q, then more power should achieve larger secrecy rate. In other words, 
we should use the maximum power. 

2) Statistical information on at the transmitter: In this case, we deal with the situation y S^- 
Taking the derivative of C S (Q) with respect to p, we get: 

9c s (Q) _ i r i 



9 P P hro ll + phl„Sg /2 QSg /2 K 



- iE h „ j — \ — — ) . (43) 

Here, we have used (T58T ) and d59l ) in Appendix |A] According to Ostrowski theorem |[24l p. 224], we 
know that if A y B and B y 0, then A fe (A 1 / 2 QA 1 / 2 ) > A fe (B 1 / 2 QB 1 / 2 ), where A fe (-) denotes the fcth 
eigenvalue arranged in decreasing order. Since ^ Eg, similarly to the methodology in Appendix lAl 
it is easy to prove that 9C $^ > 0. Thus, more power should achieve larger secrecy rate. In other words, 
we should use the maximum power. 

Remarks: For the situation ^ whether or not C S (Q) > imply that 9C q^ > has not been 
proved. This can be our future work. 

VI. The Optimal Q Under High SNR 

In this subsection, we give an analysis for high SNR, i.e., p — > oo. Our results show that for high SNR, 
the optimal Q has rank one, i.e., beamforming is optimal. The detailed analysis is given as follows. 

A. Full CSI about h/j at the transmitter 

According to E jlV-AI the optimal Q always has rank one. Let Q = uu^ with ||u|| 2 = 1. 
For high SNR, by using the fact that log(l + x) ~ log a; for large x, we write 

C a (Q) « logCphJjUuthfl) -E hB {log(p4WM} 

= logCuthflh^u) -E h „{log(htS^ /2 uutsi /2 h„)} 
= log(uth R h^u) - E ftBil {Iog(utSi;u|Vi| 2 )} 



logCuthflhJju) - log(ut££u) - E{log | Vi I 2 } 
uth/jhjjii 



-„-ln- 

lo S — + ^ + 7 



<log(h{ i E- 1 h i i) + 7 (44) 
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where we have used the fact that XJh w and h w have identical distributions for any unitary matrix U, 
and E{log l/i^i] 2 } = —7 where 7 = 0.577216- •• is the Euler's constant (since 2|/i tU) i| 2 ~ X 2 (2), 
i.e., the chi-square distribution with degree of freedom 2). In (l44l ). the maximum is achieved when 

U = V- E 1 hR/\\V E 1/2 h R \\. 

B. Statistics information about h/j at the transmitter 
For high SNR, similarly, we write 

C S (Q) 

« EhJlog^Qh^)} -E hB {log(p4Qh B )} 
= E^{log(hts|Qs|h w )}-E hw {bg(ht J s|Qs|h t0 )} 
, , a k \h w k \ } - E hm {log } b k \h W)k \ } 

k=l ' * — 'fc=l 

= E h /log S 1 ?^!! ) (45) 



11 11 

where a&'s and fr^'s are the eigenvalues of SJjQSJj and XJ,QXJ, arranged in decreasing order, respec- 
tively, and we have used the fact that Uh^, and have the identical distributions for any unitary matrix 

1 1 

U. Let pi and p nT be the largest and smallest eigenvalues of Xj^X^ Xj^. By writing 

4q4 = 4 s ? (4q4) s z4 ( 46 ) 

and applying Ostrowski theorem ll24l p. 224], we have 

a k = b k 6 k , p nT <0 k < px(k = 1,- •• ,n T ). (47) 

ii 11 

Since rank(EjjQEjj) = rank(S|,QSg), the number of non-zero elements of a k and b k s is the same, 
denoted by n. Thus, from d4"71 ), we have 

ak 



max 

k=l,-,n b k 



<Pi = A max (S|S^ 1 S|). (48) 



To proceed, we need the following lemma. 

Lemma 8: For x k > 0, y k > 0, & = !,••• , n, if /iofc?s 



2/1 H VVn k=l,-,n y k 

The proof is simple: it is based on the following 



Xl + '" + Xn < max ^. (49) 



— ^— < max {-,-), Va, 6, c, d > 0. (50) 

c + a c d J 
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From (|48T ). using Lemma [U we get that: for any h„, ^ 0, 

En T it |2 v-M! it |2 

k=l a k\iT'W,k\ _ 2-^k=l a k\ w ,k\ 

YJk=i^k\h w ,k\ 2 YJk=i h\h w> k\ 2 

, a k 
< max — 

fc=l,- - ,n Ok 

(51) 

In (1511 ). the maximum is achieved simultaneously for any h„, 7^ when Q = uouj, uo = S B 2 xo/||S £ , 3 xo| 

with xo being the eigenvector associated with the largest eigenvalue of ~E E 2 XrX!^ 2 , and correspondingly, 

ai = ujs^uo, 61 = ujs^u , 02 = • • • = a„ T = 0, 62 = • • • = b nT = 0, i.e., n = 1. To see why this is 

1 j_ _j_ _i 

the case, noting that A max (5]^5] £ ; 1 5]|j) = A max (5] £ ; 2 S/jS^ 2 ), it is easy to verify that 

f = x og^-x = A _ (s -| Ei?s -f } (52) 
01 xjx 
The desired result follows. Now, combining (l45l) and (I5TT ). we have 

a(Q)«E h llogSl^l^J! 

I Z^fc=l Ofc|«uj,fe| 

< log(A max (s|si 1 E|)). (53) 
In (l53l . the maximum is achieved when Q = uouj. Thus, the optimal Q has rank one. 

VII. Fixed Point Iteration 

In this section we propose an algorithm to search for the KKT solution according to Theorem Q] 
When and Q commute, + jI nT and Q commute for any real number 7, and vice versa. Let 
7 = (1 + 0) max{0, -A min (0)}, (3 > and let K = + jI nT . It holds that K y 0. From fl2i we get 

KQ = Tr(KQ)Q. (54) 

Equation d54l looks like the eigenvalue equation Ax = Ax, where Tr(KQ) is the eigenvalue and Q is 
the corresponding eigenvector. Recall that the power iteration method is a classical method for computing 
the eigenvector associated with the largest eigenvalue of a matrix ||2T1 p. 533] 

Xfc+i = ^ Xfc ,fc = 0,l,--- ■ (55) 

1 1 -^- X A; I 



We can derive the similar algorithm. Note that there is a difference between (1541) and the eigenvalue 
equation: Q is a Hermitian matrix. Thus, the iteration (l55l ) cannot be used directly. From (l54l . since K 
and Q commute, thus, we have that KQ = K 1 / 2 QK 1 / 2 and hence 

K 1 / 2 OK 1 / 2 * 

Q = - fm <56) 
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Note that /(Q) >z and Tr(/(Q)) = 1 for any Q € Q. The equation (l56l) defines a mapping from a 
convex set to itself: U — > O, Q i-> /(Q). The optimal Q corresponds to a fixed point of /(Q), i.e., 
/(Q°) = Q°. To search for the KKT solution, the iterative expression is 

Q k+1 = f(Q k ),k = 0,l,--- (57) 

The initial point Q° can be set to I nT /n^, or any Q° € fi. The iterations stop when the relative error 
of C S (Q) in the successive iterations is less than a preset value, e.g., 1CP 3 or 1CT 6 . If the convergent Q 
satisfies (TT3T) . we obtain a KKT solution, otherwise, we choose a different initial point. 

VIII. Numerical Simulations 

In this section we provide some examples to illustrate the theoretical findings. We assume that is 
normalized as Tr(E^) = tit, and Y^e is multiplied correspondingly by a factor -q. In simulations, we 
assume that the correlation matrices of legitimate and eavesdropper channels follow the Jakes' correlation 
model |fT8H . i.e., for p, q = 1, • • ■ , tit 

^r(p,q) = Jo{<Pr\p - q\2ird/\), 
V E (p,q) =vJo{<pE\p-q\2ird/\) 

where Jo(-) is the zero-order Bessel function of the first kind, d is the element spacing, A is the 
wavelength, and (or 4>e) is a parameter that controls the correlation among antennas and has its 
value determined by the distance between the transmitter and receiver, and the incident angle of the 
wavefront. We set d/X = 1/2. 

A. The transmitter has full information about the legitimate channel and only statistical information 
about the eavesdropper channel 

We consider a MISO wiretap channel where tit = 4, ur = tie = 1- We set h# = [0.4282 + 
0.0403i, 0.8956 + 0.6771i, 0.7310 + 0.5689i, 0.5779 - 0.2556i] T and <J) E = 0.3, rj = 0.3. The eigenvalues 
of hfihjj - S E are 2.0946, -0.0020, -0.1584, -0.4315. 

Fig. |2] depicts the function C s (z) defined in ([20]> for z in [0.01, 0.99] with step 0.01 and SNR = 10 dB. 
Among these points, the optimal point is (0.55, 2.8413) also depicted in Fig. |2] 

Fig. [3] depicts the ergodic secrecy rates during the iteration of the algorithm of Section Will for 
SNR = 10 dB and Q° = jl^. The convergent ergodic secrecy rate is 2.8413. We can see that the 
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algorithm converges rapidly. If we do 300 iterations for SNR = 10 dB and Q° = jl^, the convergent 
values are: 

A(Q) = {1.0000,0.0000,0.0000,0.0000}, 
A(0) = {0.4385, -0.0105, -1.3155, -2.2006} 
Tr(0Q) = 0.4385. 

Fig. [4] plots the ergodic secrecy rates for different 4>E from 0.2 to 0.9. It can be seen from Fig. [4] that 
the ergodic secrecy rate decreases first, and then increases with 4>e- Fig. [5] plots the ergodic secrecy rates 
for different SNR and 4>e = 0.3. As is revealed in EjV] the ergodic secrecy rate increases with SNR. 

B. The transmitter has only statistical information about both the legitimate channel and the eavesdropper 
channel 

We set <j) R = 0.5, (j) E = 0.3, t] = 0.3. The eigenvalues of Y, R - S E are 1.3503, 0.9848, 0.4432, 0.0217. 

Fig. [6] depicts the ergodic secrecy rates during the iteration for SNR = 10 dB and Q° = II4, while 
Fig. 13 depicts the ergodic secrecy rates for 30 random Q° G fl We can see that the algorithm converges 
rapidly. If we do 300 iterations for SNR = 10 dB and Q° = ^4, the convergent values are: 

A(Q) = {0.5129,0.4871,0.0000,0.0000}, 
A(0) = {0.7452, 0.7452, -2.5972, -5.4438} 
Tr(0Q) = 0.7452. 

We can see that the convergent Q has rank two. 

Fig. [8] plots the ergodic secrecy rates for different cpR from 0.4 to 0.9. It can be seen from Fig. [8] 
that the ergodic secrecy rate increases with 4>r. Fig. [9] plots the ergodic secrecy rates for different </>g 
from 0.2 to 0.4. Fig. [10] plots the ergodic secrecy rates for different SNR. As is revealed in £jV] when 
y Eg, the ergodic secrecy rate increases with SNR. 

IX. CONCLUSION 

We have investigated the problem of finding the optimal input covariance matrix that achieves ergodic 
secrecy capacity subject to a power constraint. We extend the existing result to nontrivial covariances of 
the legitimate and eavesdropper channels. We have derived the necessary conditions for the optimal input 
covariance matrix in the form of a set of equations and propose an algorithm to solve the equations. 
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Appendix A 
Proof of Lemma Q] 

We prove the result in two parts. First , we prove that if E# — E# y 0, then C s < 0. Since h# ~ 
C7V(0, T, R ), h E ~ CM(0, E B ), we can write 

h/j = S^X, (58) 
h £ = E^X- (59) 

By inserting (I58T ) and (|59l into ©, we get 

C S (Q) = E h „{log(l + pht^QEj/X)} 

- E h Jlog(l + ph^E^E^X)}. (60) 

1/2 1/2 1/2 1/2 

Let xi > x 2 > ■ ■ ■ > x nT and yi > y 2 > • • • > 2M T be eigenvalues of E^' QS R ' and E^ QE^ , 
respectively. By using the fact that Uh^ and have the identical distributions for any unitary matrix 
U, we have 



C S (Q) = E h J log (l+p£ ; 

^ ^ i=i 

-E hit ,|log^l + p^ yi |Vi| 2 )}- (61) 

According to Ostrowski theorem (2H p. 224], we know that if A >z B and B y 0, then A fc (A 1 / 2 QA 1 / 2 ) > 
A,t(B 1 / 2 QB 1 / 2 ), where A^(-) denotes the kth eigenvalue arranged in decreasing order. With this, we know 

that Xi < yi, i = 1, • • ■ , tit- On the other hand, it is easy to verify that the following function 

c / nT 

g(z 1 , ■ ■ ■ ,z nT ) = E 



h„|log + (62) 



is strictly increasing with respect to Zi, i = 1, ■ • ■ ,ny. Thus, we get that C S (Q) < for any Q. This 
completes the first part. 

Second , we prove that if Er — Eg is none negative semi-definite, then there exists a Q such that 
C S (Q) > 0. Let u be the eigenvector associated with the largest eigenvalue A of E# — Eg. Since 
A > 0, we get that uT(E# — Se)u = A > and u^Eru > u^S^u. We will prove that Q = uu^ 
achieves C S (Q) > 0. In this case, we know that x\ = u^E^u, x% = ■ ■ ■ = x nT = 0, y% = u^E^u and 
y2 = ■ ■ ■ = Un T = 0. Since x\ > y\ and the function g{z\, ■ ■ ■ , z nT ) defined in d62b is strictly increasing 
with respect to z\, we get that C S (Q) > 0. This completes the proof. 
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Appendix B 
Proof of Lemma [2] 

From ® and d59l ), we get 

C s (Q) = log(l + ph^Qh R ) 

- E hu ,{log(l + ph+^QE^X)}- (63) 

It follows from the Jensen's inequality [26, p. 25] that logE(x) > E(logx). With this and the fact that 
EhJh^Ah^} = E h jTr(Ah,X)} = Tr(A), we get 

C S (Q) > log(l + ph^Qhn) - log(l + P Tv(Q^ E )). (64) 

Note that hJjQh/j— Tr(Q5]^) = Tr(Q(h^hJj— £#)). Let u be the eigenvector associated with the largest 
eigenvalue A of hRh^ — Sis. Since A > 0, we know that Tr(uu^(h^h|j — ~Ee)) = ut(h/jh^ — Xl^u = 
A > 0. Thus, C S (Q) > holds for Q = uut This completes the proof. 

Appendix C 
Proof of Lemma [3] 

From |[TT1 Eq. (64)], we know 

E z {log(l + pztQz)} 

00 er l ( 1 A 

1 " ~i TTT - '. TTTTtT^TTToT dt. (65) 



E z {log(l + pztQz)}= / e -*(---_ 3r ___)dt. (66) 



/o t V det(I + tpRV2QRi/2 
By inserting R 1 / 2 QR 1 / 2 = UiDiU| into (g5]l, we get 

"OO /J 

/o \J + 
Performing partial fraction expansion, i.e., 

1 1 Pdj e 

* tn£i(i+fA*) ^nS,-(i-*/4f) 1+t ^ 

and using 

oo -t y 

——dt=-F 1 (a), (68) 
1 + ta a 



o 



we get ©. This completes the proof. 
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Appendix D 
Proof of Theorem Q] 

It follows from (© that \I/Q = Q\I/ = 0, that is, \I/ and Q commute and have the same eigenvectors 
p. 239] and their eigenvalue patterns are complementary in the sense that if Aj(Q) > 0, then Aj(\l/) = 0, 
and vice versa lfl4l . This result, when combined with ([8]), implies that and Q commute and have the 
same eigenvectors. Further, we get 0Q = Q0 = 9Q, which, when combined with Tr(Q) = 1 and the 
fact Tr(Q0) is always real, leads to 6 = Tr(Q0) and £[2]) (also see ifBTl ). 

The condition (ITZb reveals that for the optimal Q, Q0 is a scaled version of Q. Further, the eigenvalues 
of corresponding to the positive eigenvalues of Q are all equal to Tr(Q0), while the remaining 
eigenvalues of are all less than or equal to Tr(Q0), which follows from ([8]), (fT2l) and \I/ >z 0. Based 
on the above (fT3l) follows. 



Appendix E 
Proof of Lemma H] 

Denote the expectation in the left hand side of (fl4l as /. We write z = ~R}/ 2 h w where h n 
[h Wt i, ■ ■ ■ , h WjnT ] T ~ CJ\f(0,l nT ), and /i^'s follow i.i.d. CM(0, 1). With this, we have 

/ 9R 1 /2 h „,htR 1 /2 



I = E h J r- -u,-^- , (69) 

ll + phtRl/2QRl/2 hw J 

By inserting R^^R 1 / 2 = UiDiUj into d69]), we have 

/ = P R 1 / 2 E h ( h w ht | Rl/2 

1 1 + phLUiDiUih^ J 

Then we use the fact that Uh^, and h w have the identical distributions for any unitary matrix U to obtain 

/ = pR 1/2 UiYUlR 1/2 (71) 



where 



with (i,j)th entries given by 



h h T 
ll w ll w 



Y = E h J ^ (72) 

1 + phLDxh, 



1 + pJ2k=i d k\K>,k 



From the gamma integral [16] we have 

1 1 f°° 

— = ^tt / t z - l e~ ta &t, Re(js) > 0,o > (74) 
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where T(z) = J °° u z x e u du, we let z = 1 to obtain ^ = J °° e * a di. With this identity, we can write 



r°° ( . . 1 

Y ij = ^ e- t E hw |fc w , < ^JJe-^*l fc -*l a |dt. (75) 

Since h Wt k's follow i.i.d. CA/"(0, 1), we know Yjj = for i ^ j, i.e., Y is a diagonal matrix with (k, k)th 
entries given by 

f°° e~* 1 
Ykk = / TTaTTi 7TT— 7TT dt ' k ^ M (76) 



Y fefc = / ^j-: __dt,A;>M. (77) 



/o n-=i(i+^ 

These integrals can be easily calculated. Performing partial fraction expansion 



i A i i 



_y- L 

n£x(i + « ^ nfi-(i - 1 + ^ 

and using (l68l ) and 



(78) 



00 ^ d * = m^im (79) 

o (l + ta)(l + t&) a-b K ' 

2 dt = F 2 (a), (80) 



oo e -t 







;i + ta) 2 



we get (fl~4l) . This completes the proof. 



Appendix F 
Proof of Lemma [5] 

We need to prove that 

4>{tz x + (1 - t)z 2 ) < t^zx) + (1 - t)(j)(z2), Vt G [0, 1]. (81) 

Let Xi and X2 be the optimal X associated with Z\ and z%. Consider the problem associated with 
tzi + (1 - t)z 2 , i.e., 

min Tr(££X) (82) 
s.t. Tr(h R h^X) = (tz! + (1 - t)z 2 )\\h R \\ 2 , 
Tr(X) = 1, 

x y 0. 

It is easy to verify that tX.\ + (1 — t)X. 2 satisfies the constraints in the problem of d82l ) with the 
corresponding objective value t(j)(z\) + (1 — t)cf)(z 2 ). Thus, ( f8TT) holds and 4>{z) is a convex function. 

March 12, 2010 DRAFT 



20 



Appendix G 
Proof of Lemma [6] 

Let A > be an eigenvalue of aa^ — A, and we have 

det(aa f - A - AI) = 0. (83) 

Note that A + AI is positive definite. By using the identity det(B — aa^) = (1 — a^B *a) det(B) for a 
positive definite matrix B, it follows from (l83l that 

1 - a f (A + AI) _1 a = 0. (84) 

Denote £(X) = 1 — a"!" (A + AI) _1 a. It is easy to verify that £(X) is a strictly increasing function. Thus, 
£(A) has only one positive root, and ^(0) ^ 0, i.e., is not a eigenvalue of aa^ — A. Thus, all other 
eigenvalues are negative. This completes the proof. 

Appendix H 
Proof of Lemma [7] 

We have 

log x > E(log Y) => E (log > 0. (85) 



By using Jensen's inequality 11261 p. 25], we have 

logE(^) >E(log^ 

=>;<e(p). <8« 
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Fig. 5. Ergodic secrecy rate for different SNR. Full CSI on h_R is used. 
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Fig. 6. Ergodic secrecy rate during the iteration. Only statistical info, on hft is used; Q° = jLi, SNR = 10 dB. 
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Fig. 7. Ergodic secrecy rate during the iteration. Only statistical info, on hft is used; 30 random Q°, SNR = 10 dB. 
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